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ABSTRACT. We provide an axiomatic framework for Quantum Field Theory at finite 
temperature which implies the existence of general analyticity properties of the n-point 
functions; the latter parallel the properties derived from the usual Wightman axioms in the 
vacuum representation of Quantum Field Theory. Complete results are given for the prop- 
agators, including a generalization of the Kallen-Lehmann representation. Some known 
examples of "hard-thermal-loop calculations" and the representation of "quasiparticles" 
Q\ . are discussed in this general framework. 



D . 1. INTRODUCTION 



In the past fifteen years, there has been an increasing interest for the physics of media 
at very high temperature; in the latter, the basic quantum fields of matter are supposed 
to manifest themselves through new properties which highly deserve to be investigated 
from a theoretical viewpoint. In particular, the so-called "quark-gluon plasma" has been 
the object of a number of theoretical studies and computations [1-3]. These are generally 
based on an adaptation of the standard methods of quantum statistical mechanics to the 
relativistic system consisting of the basic fields of quantum chromodynamics. The fields 
are supposed to be in an equilibrium state of infinite volume at temperature T = 1/(3 in a 
certain privileged Lorentz frame; the latter fixes the time and space variables (xq, x) and 
the corresponding Fourier conjugate variables, namely the energy and momentum (uj,p). 

The use of the Matsubara imaginary-time formalism (ITF) (see e.g. [4] and references 
therein) has resulted in computations involving discrete imaginary-energy summations for 
obtaining perturbative approximations of the retarded propagators of the fields in the space 
of (complex) energy and (real) momentum variables. In this formalism, the quasiparticles 
are associated with "modes" obeying a (real or complex) "dispersion law" uj = f(p), which 
appear as poles of the form Z{j>)/[ui — f(p)] in the retarded propagators of the fields. Such 
a structure has been displayed in the so-called "hard-thermal-loop calculations" which 
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aim to exhibit the very high temperature behaviour of a quark-gluon plasma; summing 
the one-loop self-energy contributions in a leading approximation at large T yields the 
following result: the gluon propagator exhibits two typical poles uj = ft(p) and uj = fi(p), 
interpreted respectively as "transverse and longitudinal plasmon modes" [5,6]. 

The independent elaboration of a formalism which only involves real-time quantities 
(RTF) has resulted in a double-field matrix formulation, now called thermo-field dynamics 
(see e.g. [7] and references therein). This approach introduces as preferable basic functions 
the time-ordered and anti-time-ordered expectation values. 

In the past years, there has been a long debate about the consistency of the ITF and 
RTF approaches, the latter being both expressed in terms of appropriate versions of the 
path-integral formalism; the controversy has generally been set in this framework, where 
genuine subtleties appear in handling "time-ordered paths" in the complex plane of the 
time- variable (see e.g. [8]). Although the situation now seems to have been clarified in 
favour of the consistency [9], one feels a real need for a general structural study of the 
thermal Green functions going beyond the time-ordered path technique. 

In the present work, we adopt such a general (non-perturbative) viewpoint on Ther- 
mal Quantum Field Theory: we suggest that the validity and consistency of the ITF and 
RTF formalisms are by-products of model-independent structural properties of the ther- 
mal n-point functions which follow rigorously from an appropriate set of general physical 
principles. As a matter of fact, we will show that one can define an axiomatic program 
for studying the analytic and algebraic properties of the n-point functions of fields in 
a thermal equilibrium state. This program is completely similar to the one which has 
been developped in the sixties for Quantum Field Theory in the vacuum state within the 
Wightman axiomatic framework and which has led to such important structural insights 
as the PCT and Spin-Statistics Theorems, collision theory, dispersion relations and the 
Osterwalder-Schrader Theorem (see [10-13] and references therein). 

The main discrepancy of our axiomatic starting point with respect to the familiar 
Wightman axioms will consist in releasing the Lorentz covariance properties and replacing 
the spectral condition by an appropriate formulation of the KMS- condition: this idea relies 
on the basic analysis of [14], completed by the results of [15,16] which display the KMS- 
property as being a general criterion for systems in equilibrium states. Moreover, in spite 
of the breaking of Lorentz covariance due to the thermal bath, the basic role played by 
the causality cone V + = {x = (xo,x);xq > \x\} for quantum field systems implies that a 
(stronger) relativistic form of the KMS-condition can be justified [17], as a remnant of the 
relativistic spectral condition. 

In this preliminary work, we shall present substantial results for the structure of the 
two-point functions; only partial results will be given as far as the general program of 
n-point functions is concerned. However, the case of the thermal two-point functions 
being already of physical interest, we give a special importance to a Kallen-Lehmann-type 
representation (already presented in [18,19]) which we are able to derive from the general 
principles; among other advantages, this representation opens a new possibility for the 
characterization of particles in Thermal Quantum Field Theory. 
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After having presented our axiomatic framework for the representations of Quantum 
Field Theory in thermal equilibrium states in Sec. 2, we devote Sec. 3 and 4 to the pure 
implications of locality: this part of the program reproduces in the setting of Thermal 
Quantum Field Theory properties of general field theory in the vacuum state. It concerns: 

a) in Sec. 3, the basic analyticity properties of n-point Green functions in the space of 
complex energy and momentum variables, 

b) in Sec. 4, the derivation of Kallen-Lehmann-type integral representations for the two- 
point commutator and retarded functions which are valid in the absence of Lorentz covari- 
ance and of spectral condition. 

The implications of the KMS-condition will be studied in Sec. 5. In particular, we 
shall exhibit the resulting double analytic structure of the two-point functions in the time 
and energy variables and the relations between the Fourier transforms of the retarded 
and time-ordered functions which replace the usual "coincidence relations" implied by 
the spectral condition. The derivation of similar properties for the n-point functions will 
only be initiated there and mentioned as a further important part of our program to be 
implemented. 

In Sec.6, the following complements on the structure of two-point functions will be 
presented: 

a) incorporating the KMS-condition into the Kallen-Lehmann-type formula of Sec. 4 in 
order to provide a corresponding representation for the two-point correlation function 
itself, 

b) completing the results of Sec. 5 by analyticity properties in the complex space and time 
variables which result from our relativistic form of the KMS-condition, 

c) studying the basic Feynman-type operations on two-point functions and illustrating 
all the previous results on the perturbative examples mentioned at the beginning of this 
introduction, 

d) giving a comparative discussion of two different characterizations of the notion of particle 
in Thermal Quantum Field Theory. 



2. THE AXIOMATIC FRAMEWORK 

The adaptation of the Wightman axiomatic framework to the case of quantum fields in a 
thermal equilibrium state Vtp of temperature T = 1/(3 relies on the following well-known 
ideas. For simplicity, we consider the case of a single hermitian field (j)(x), namely a sys- 
tem of "observables" <f>(f) = f <p(x)f(x)dx depending continuously on test-functions / on 
Minkowski space (taken for convenience in the Schwartz space <S(1R 4 )). These observables 
form an algebra which satisfies the general axiom of "locality" or "local commutativity" , 
namely: [</>(/), <f>(g)] = for all pairs of test-functions (f,g) whose supports are mutually 
space-like separated in M 4 . 

This axiom, which expresses the principle of Einstein causality, is independent of the 
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representation in which a system is described; it therefore holds in the Hilbert space Hp 
of thermal states in which the field observables </>(/) are supposed to act, as operators 
defined on a suitable dense domain generated by the state vector Qp. It is understood that 
Qp determines a fixed Lorentz frame, i.e. a distinguished set of time and space variables; 
the corresponding unitary representation Ut of the time-translation group in Tip (i.e. the 
evolution operator group) leaves Qp invariant, since an equilibrium state is stationary. We 
shall only consider here the case when Qp is also invariant under the (unitary) representa- 
tion of space-translations Us in Tip. The action of the space and time-translations on the 
field is, as usual: 4>{xq + t, x + a) = UtUs<f>(xo, x)U^ 1 U^ . The previous axioms do not 
differ from the corresponding Wightman axioms in the Hilbert space H vac generated by 
the vacuum state fl vac ', however, there are no unitary operators which implement Lorentz 
transformations in Hp. 

In order to complete our axiomatic framework in Hp, we shall introduce the "corre- 
lation functions" or "Wightman functions" Wn(xi, . . . ,x n ) =< Qp, 4>{x\) • • -cf)(x n )Qp > 
and adopt the viewpoint of the reconstruction theorem [10]. A Thermal Quantum Field 
Theory is thus entirely specified by the knowledge of the corresponding set of tempered 
distributions {Wi/ 3 "*; n G N} which have to satisfy the usual positivity conditions of Wight- 
man functions [10]. (It can equivalently be presented as the representation associated 
with a certain positive functional on the Borchers-Uhlmann algebra [20,21] of terminating 
sequences of test-functions {fo, fi(xi), • • • , fn(xi, • • • , x n ), 0, . . .}.) 

Apart from the properties which express the previous axioms, namely: 

i) Locality: 

(xi, . . . , Xi, x i+ i, ...,x n ) = (xi, . . . , x i+ i, Xi...,x n ) (1) 

for Xi - x t+1 £ V+ U W with V- = -V+, 

ii) Translation Invariance: 

WP(x 1 ,x 2 ,...,x n ) =WP(x 1 + a,x 2 + a,...,x n + a) (2) 

for all a in IR 4 , the thermal Wightman functions should satisfy the following analyticity 
properties which replace those implied by the spectral condition in the vacuum case. These 
properties are the mathematical expression of the fact that the state Qp is in thermal 
equilibrium [14-16]. 

iii) KMS-condition: 

Let eo be the unit vector of the time-axis. For every n G N and every pair (/, J) of 
ordered sets I = {l,...,m}, J = {m + l,...,n}, there exists a tempered distribution 
Fu(xi, . . . , x n ; t) which is holomorphic with respect to the complex variable t in the strip 
—(3 < Imt < and admits distribution boundary values on {t G R} and {t G R — i(3} (still 
denoted by Fjj) satisfying the following conditions: 

Fjj(xi, . . . , x n ; t) = ({x t + te } ie i, {xj}j€j), (3) 
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F IJ (x 1 , . . . , x n ; t - iff) = W$?\{xj} je j, {xi + te } ieI ). 



(4) 



Let us introduce the Fourier transforms of the Wightman functions, which (in view 

of ii)) can be written as W$f\pi, ■ ■ ■ ,p n ) 5(pi H \-Pn), with pi = (u>i,pi),l < I < n. 

Then the distributions Wn^ on the linear manifold p± + ■ — h p n = satisfy the following 
property which is equivalent to iii): 

iii') KMS-condition in the energy variable: 

For each pair (I, J), the following identity holds: 

wW(J,I) = e-^'WW{I,J), (5) 

where u r = J2iei u * and ^^(I, J) stands for ({Mie/, {pj} je j)- 
Remark: If we introduce the commutator functions 

C^\I,J)=W^(I,J)-W^(J,I) (6) 

it follows from Eq.(5) that VV^ is "essentially" determined from (i.e. determined up 
to a part which factor izes S(ui)) by the formula: 

(l-e-^)W^(I,J)=C^(I,J). (7) 

We notice that this formula replaces the usual one expressing the spectral condition in the 
case of the vacuum state: 

W^ ac \l,J) = 0{oo I )C^ c \l,J) (8) 
(9 denoting the Heaviside step-function). 

In the limit of zero temperature {(3 — > oo), Eq.(7) tends to Eq.(8); the "inverse Bose- 
Einstein factor" (1 — e _/3w/ ) in Eq.(7) implies (since VVn^ is a tempered distribution) that 
Wn(I, J) has an exponentially decreasing behaviour of the form e - ^'^ 7 ' for ui tending to 
— oo, instead of being equal to zero for ui < as VVn a °\l, J)- The occurence of negative 
energy contributions in corresponds to the possibility of "extracting energy from the 
ambient thermal bath" . The exponential suppression of these contributions is a remnant 
of the condition of positivity of energy in the vacuum case. In fact, in the relativistic 
theory the latter is true in all Lorentz frames (relativistic spectrum condition); this implies 
a more stringent support condition for Wn ac \l, J) than the condition ioj > exhibited 
by Eq.(8), namely: 

supp Wi vac \l, J) C {(pi, . . ■,p n )',Pi = -P. j e V+}, (9) 
where pi = J^ieiPi and P J ^ s defined analogously. 
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A corresponding relativistic form of the KMS-condition has been justified in the anal- 
ysis of [17] based on the general principles of "Local Quantum Physics" [22]; this analysis 
gives a firm background to the following properties (stronger than iii), hi')) of the Wight- 
man functions W$?\l, J). 

iv) Relativistic KMS-condition: 

For each pair (I, J) , the distribution Fjj of condition iii) admits an analytic continuation 
(which we still call) Fu(xi, . . . , x n ; z) with respect to the complex four-vector variable 
z = (zq,z) in the following tube-domain: Tp = {z G C 4 ; Imz G V~ , Imz + f3e G V + }. 

The boundary value equations (3) and (4) can then be replaced respectively by: 

lim F I j(x 1 ,...,x n ;x-ie)=W^({xi + x} ieI ;{x j } je j), (10) 

lim F IJ (x 1 ,...,x n ;x-iP + ie) = WP({x j } j€J ;{x i + x} ieI ). (11) 

In the space of energy and momentum variables, condition iv) can be equivalently expressed 
by iii') supplemented by the following condition iv'): 

iv') Essential support conditions in p-space: 

For each pair (I, J), the distribution VVn(I, J) admits the cone {(pi, ■ ■ ■ ,p n )',Pi = ~Pj £ 
V + } as a majorant of its "essential support in the sense of exponential decrease". A more 
precise formulation of this condition is that the products of the tempered distribution 

W^\l, J) by the (smooth) functions e ¥ ~ UJi ) anc [ e _ ^ W/ are required to be also 

tempered. 

For completeness, we should also mention an additional postulate whose role is to 
restrict our framework to the case of equilibrium states Vtp which are pure phases. 

v) Time-clustering postulate: 

For each pair (I, J), one has in the limit \t\ — > oo (with the notation used in Eq.(3)): 

F IJ (x 1 ,...,x n ;t) = W^({x i } ieI ) wf_ ) m ({x J } j6J )+o(l). (12) 



3. ANALYTICITY PROPERTIES 
IN THE ENERGY AND MOMENTUM VARIABLES 

Recently, a procedure for studying the usual n-point retarded and advanced functions and 
possibly introducing generalized retarded functions characterized by various analyticity 
properties in the energy variables has been presented in the traditional imaginary-time 
formalism of Thermal Quantum Field Theory [23] . Here, we would like to emphasize that 
a general understanding of the algebraic and analytic properties of these Green functions 
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follows from the Wightman axiomatic approach of quantum field theory, and that the 
latter allows one to control exactly which results are identical to those of the vacuum case 
and which ones are different. 

We devote this section to the introduction of n-point thermal Green functions enjoy- 
ing analyticity properties in the energy and momentum variables. Their definitions and 
algebraic study in the previous axiomatic framework and the (correlated) derivation of 
their analyticity properties as pure consequences of the axiom of locality are completely 
identical to those which have been given for the vacuum representation of Quantum Field 
Theory. This structure has been the object of a number of works in the past [24-28] . 

The basic idea is that for each n, one can define certain privileged combinations 
r a of permuted n-point Wightman functions multiplied by appropriate products of step- 
functions of the (differences of) time-variables. Each distribution r a enjoys remarkable 
support properties which are implied by locality, namely the support of r a is contained in a 
(salient) convex Lorentz invariant cone T a of the space of vector differences Xi — Xj. Then, 
in view of a basic result of complex analysis (see e.g. [10,29]), the Fourier transform of r a 
is the boundary value on the reals of a holomorphic function f a in a tube-domain T a of the 
complex energy-momentum space which is of the following form: {k = (k\, . . . , k n ); k\ + 

h k n = 0; Irak G r a }, where r a (called the basis of the tube T a ) is the dual cone of r a 

(i.e. the set of vectors q such that qx = q\X\ + ■ — h q n x n is positive for all x in r a ). Each 
of these distributions r Q (or holomorphic functions f a ) will be called a generalized retarded 
function. It is the expectation value < Qp, R a Qp > of a corresponding generalized retarded 
operator R a (x\, . . . , x n ). The generalized retarded operators (resp. functions) satisfy two 
sets of basic algebraic relations, namely: 

a) discontinuity relations between the various R a (resp. r a ); in this connection, generalized 
absorptive parts are introduced; 

b) relations between each R a and the (anti) time-ordered operator products. 
The two-point function: 

Let C(x) be the commutator function, i.e., C{x\ — x<i) = W?, £2) ~ W2 (#2, x\). 
Locality implies that supp C C V + U V~. There are two distributions r a in this case, 
namely the "retarded and advanced functions", defined formally as r(x) = i9(xq)C(x) and 
a(x) = —i0(—xo)C(x). 

The relation r — a = iC corresponds to the splitting of the support of C into its two 
convex components, since supp r C V+ and supp a C V~. (Note that this splitting is 
defined up to a distribution with support at the origin). 

This splitting is in turn equivalent to the following one for the Fourier-transformed 
quantities: iC(p) = r(p) — a(p), where C is usually called spectral function [1-8] and f and 
a are the boundary values of holomorphic functions (denoted similarly by) r(k) and a(k) 
(k = p + iq) in the respective tube-domains T + = {k = p + iq; q G V + } and T~ = — T + 
of C 4 . Conversely, given a holomorphic function in T+ U T~ (with at most a tempered 
behaviour at infinity and near the reals), the spectral function obtained by taking the 
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difference of the corresponding boundary values defines a commutator function which does 
satisfy locality. According to the common use, the holomorphic function r(k) (or a(k)) 
is then called the propagator whose associated spectral function is C(p). (In view of the 
Hermitian character of the field, C(p) is also the imaginary part of r(p) and is commonly 
described as such, although its characterization as the discontinuity of the holomorphic 
function f(k) across the reals is more substantial). 

Let r(x), x = x\ — X2 (resp. t(x)) be the two-point time-ordered (resp. anti-time- 
ordered) function. The following relations hold: 



r = i( T - W) = -i(r - W), a = i(r - W) = -iff - W), 



(13) 



where W{x 1 - x 2 ) =W^ ) {x 1 ,x 2 ) and W{x 1 -x 2 ) = W { 2 P > \x 2 ,x 1 ). 
Examples: 

1) The free field propagators: 

For any free scalar field of mass m, the commutator function £( m ) and therefore the associ- 
ated retarded and advanced functions A m ' and a( m ) are independent of the representation 
generated by the thermal state flp. In fact, in this case C^ m \x) is a structural function of 
the field algebra determined by the c-number commutation relations of the field. There- 
fore, for every thermal representation with temperature one has (as in the vacuum 
representation): C^ m \x) = / K 4 e~ tpx C( m Xp) dp, with 



2n 



(14) 



and 



f {m \k) (resp. d [m) (k)) 



1 



4tt< 



{kl-k 2 -m 2 ) 



(15) 



for k = (ifeo, k)=p + iqeT+ (resp. ) keT~). 
2) The Weldon-Pisarski (WP) propagators: 

The following expressions of the high-temperature transverse and longitudinal gluon prop- 
agators have been obtained by the resummation of dominant contributions at small (3 from 
the one-loop self-energy diagrams: 



A* (ft) 



kr\ k — ~ — 

° 2k 2 




(16) 



Mk) = 



k' 2 



k 2 



k + M 1 — 




(17) 



In these expressions, M denotes the "Debye screening mass" which is directly related to 
the plasma frequency. The propagators A t and A; can be checked to be holomorphic in 
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the tubes T + and T~ in the privileged sheet of the logarithm which defines the physical 
sheet of these functions. Therefore the associated two-point functions satisfy locality (in 
this connection, see [30] for a detailed numerical analysis of these examples and further 
comments on the construction of causal propagators). 

The n-point functions: 

For each n, the distributions r a are labelled by a "cell-function" a which represents a 
geometrical cell 7 a (namely a polyhedral cone) in the space = {h= (h u ...,h n ) G 
M n ; hi + - • - + h n = 0} deprived from all the hyperplanes Hj U J with equation hi = —hj = 0. 
Here, (I, J) denotes an arbitrary partition of {l,...,n} and we use the notation hi = 
^2 i&I hi. A cell a is defined as the "sign-valued function" I — > a(I) = + or — such 
that 7 a = {h G H^;WI C {1,2, ... ,n},a(I)hi > 0}. Two cells a\ and a 2 are called 
"adjacent cells along the face i/j™ j " if they only differ on the complementary sets I 
and Jo, namely if «i(/o) = — &i(Jo) = — ^2(^0) = ^2(^0)- The corresponding generalized 
retarded operators or functions will also be said to be adjacent. 

Each generalized retarded function f a is holomorphic in the tube T a whose conical 
basis r a in {q = Irak = (qi, ... , q n ); q± + ■ ■ ■ + q n = 0} is defined by the set of conditions 
a(I)qi G V + (for all proper subsets I of {1,2, ... ,n}) [27,28]. Each ordinary retarded 
function r^ 1 ' corresponds to the cell a such that ot({j}) = + for all j • E {1, 2, . . . , n} except 
for j = i. 

Discontinuity relations: 

They are generated by the following set of relations which connect all the pairs of adjacent 
generalized retarded operators [25,26]: for every partition (/, J) of {1, . . . , n} and for every 
pair of cells (cti, a^) which are adjacent along the face Hj U j, there holds: 

^Ql — Ra 2 = 2[-R a (-0, R a (J)], (18) 

where or- 1 ', denote the common restrictions of the cell- functions 011,012 respectively 
to the proper subsets of I and of J, and R a (i),R a (j) are the corresponding generalized 
retarded operators in the space of the variables labelled by the elements of I and J. 
Correspondingly, there holds the following set of relations, which generalize the relation 
r — a = iC of the two-point function: 

r ai -r a2 =i ((ftp, R a (i)R a( j)Q.p) - (ftp, R^R^Qp)) . (19) 

In the latter, the terms at the r.h.s. are interpreted as generalized absorptive parts in the 
respective channels (I, J) and (J, /). 

Steinmann relations [24]: 

The generalized retarded operators R a are not linearly independent. In fact, if two adjacent 
pairs («i, 012) and (a[, a' 2 ) along the same face Hj n J admit the same restrictions a( J \ 
the following corresponding relation holds: R ai — R a2 = R a ' — R a ^ . 
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Relations with the time-ordered operator products: 

Let T(I) denote the time-ordered product of the nj field-operators <f>(xi) for i £ I C 
{1,2,..., n}. For each cell a the following expression of R a is valid [28]: 

(-i)"" 1 R«(xi, ...,x n ) = T({1, 2, . . . , n}) + ^ (a) (-l) r_1 T(h) ■ . .T(J P ), (20) 

where the sum runs over all ordered partitions (Ji, . . . , I r ), 2 < r < n, of {1, 2, . . . , n} 
such that ct(ii) = ol{I\ U ^) = ol{I\ U . . . U I r -i) = — ■ Similar relations involving the 
anti-time-ordered operator products can also be written. Eq.(20) generalizes the relations 
(13) of the two-point function. 



4. INTEGRAL REPRESENTATION 
OF THE TWO-POINT SPECTRAL FUNCTION 

We present here an integral representation of the thermal two-point commutator and spec- 
tral functions C(x) and C(p) which is a pure consequence of locality; from a technical 
viewpoint, it also relies on the fact that C(p)/uj has to be a (positive and even) measure 
as a by-product of the positivity and KMS conditions (see Sec. 5). Our representation 
is characterized by a certain "weight-function" which plays the same role as that of the 
Kallen-Lehmann representation of the vacuum two-point function. This weight-function 
can also be reconstructed from C(x) by a simple inversion formula. 

Proposition: 

i) The following integral representation holds: 

/•OO 

C(x)= dm D(x;m) C {m \x); (21) 
Jo 

in the latter, the weight-function D(x; m) is a tempered distribution with support in M 3 xR+ 
which is uniquely defined and computable in terms of C by the following inversion formula 

D(x;m) = 2nr— (9(m) J dx x Jo(mVx^) C(x)), (22) 

where Jo is the zeroth-order Bessel function of the first kind. 

ii) The following corresponding representation applies to the spectral function: 

C(oj,p) = 1 * e(u) I du [ ds 5{u 2 - {p-uf - s) p{u,s), (23) 
(2tv)2 J r s Jo 

where 

P% *) = * 5 r I dx e-**D& y/a) (24) 
2 (27r) 2 y/s Jr3 
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belongs to S'(R 4 ), with supp pcK 3 xK+. 

Similar results were first stated by Gervais and Yndurain in some unpublished paper 
[31]; Eq.(23) can be interpreted as a Jost-Lehmann-Dyson representation (see [12] and 
references therein) in the absence of spectral condition. The proof which we present here 
is self-contained; the missing technical details of distribution theory will be given elsewhere 
[32]. 

a) Gauss-type transforms in the time and energy-variables: we first associate with the 
(4d-) Fourier pair (C,C) the following (3d-) Fourier pair (^(x; A), \&(p; A)) 

ty(x; A) = % / dxo e 4A xq C(xq, x), (25) 

*(p;A) = (2A)* J dw e~ Xuj2 u C{u;,p). (26) 

It can be checked that \1/ and \1/ are tempered distributions (resp. in x and p) which are 
(slowly-increasing) holomorphic functions of A in the complex half-plane C + = {A; ReX > 
0}. Moreover, the mapping C — > ^ and (thereby) C — > ^ are one-to-one since the integral 
transformation (26) can be considered as a Laplace-transformation in the variable u 2 and 

therefore inverted (here, one makes use of the fact that C ^\f^ is an even measure in a;). 

b) The transforms $ and We define 

$(£;A) = e£ tf(x;A). (27) 

Locality implies that the distribution- valued holomorphic function $(x, A) is tempered like 
A). By inverting Eq.(27), we obtain the following relation between the Fourier trans- 
forms $(p; A) and ^(p; A) of $ and 

3 

V(p;\)=(-Y f due- x (P- a)2 $(u;\), (28) 

V 71 "/ 

the latter being valid for all A in C + . 

c) The "weight-functions" p and p: we call respectively p(x; s) and p(p; s) the inverse 
Laplace-transforms of $ and $ with respect to the variable A; (p, p) is a (3d-) Fourier pair 
of tempered distributions on M 4 with support contained in M 3 x R+. For our needs, we 
write explicitly the mappings $ — > p and p — > 

i r +QO 

p(x;s) = — du e lus $(x;iv), (29) 

/•+oo 

$(p;A)= / dse- Xs p(p;s). (30) 
./o 
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d) The integral representations (21), (23): by plugging Eq.(30) into Eq.(28), we obtain 

3 

#(#A)=(-Y / du r ds e- x (^- a)2+s )p(u;s) 
V 71 / Jr'-i Jo 



A^ 

7T 



e- Xuj du 2 ( du ds 5(u 2 -(p-u) 2 -s)p(u;s)\ . (31) 



o 



By comparing Eqs. (26) and (31), we see that (2tt) ^9(oo)C(u,p) and the bracket in the 
r.h.s. of Eq.(31) admit the same Laplace-transform with respect to the variable uj 2 and are 
therefore equal in view of a). Since C(u,p) is odd in the variable u, Eq.(23) is therefore 
established and Eq.(21) immediately follows by Fourier transformation. 

e) The inversion formula (22): by plugging Eq.(25) into Eq.(27) and the latter into Eq.(29), 
we readily obtain (after an admissible interchange of integration, the following formula 
being understood in the sense of distibutions in the variables x and s): 

/+oo / • f+OO 2 \ 

dx x C(x) (-J du e iv8 ~^\ . (32) 



Since the bracket in the r.h.s. of Eq.(32) is an integral representation of i-^ (9(s)Jq(Vsx 2 )) 
(see [33], page 357, formula (11)), Eq.(32) takes the form (22) if one replaces p by the 
distribution D(x; m) = Airni p(x; m 2 ) (linked to p by Eq.(24)). 

From the previous proposition, one easily deduces the following 

Corollary: 

The (Fourier pair of) retarded thermal propagators r(x) and f(k) satisfy the following 
general integral representations 

POO 

r(x)= I dmD(x;m) r {m) (x) (33) 
Jo 



oc 



r (^) = ~T I du I ds =; p(u; s) , (34) 

(27f)2 y K 3 J (k* — (k — u) 2 — s) 

where k varies in T + . The same representation holds for d{k) with k varying in T~ . 

Remark: 

The Kallen-Lehmann representation, which is valid for a general propagator in the vacuum 
state in the axiomatic Wightman framework (see [12] and references therein), appears as a 
special case of the representation (21); it is obtained when the weight-function D is of the 
form D(x;m) = 2mpo(m 2 ) (i.e. p(u;s) = (2%)^ po(s) S(u)), p being a tempered (positive) 
measure with support contained in M+. 
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5. CONSEQUENCES OF THE KMS-CONDITION 

The absence of the spectral condition has a drastic effect on the structure of the n-point 
Green functions in the space of energy and momentum variables: as we shall explain it 
below, the various generalized retarded functions in complex energy-momentum space do 
not admit in general a common analytic continuation; in contrast with the case of the 
vacuum representation, there do not exist analytic n-point Green functions in (connected) 
primitive domains for the thermal representations of quantum Geld theory. It also appears 
as a related fact that the (anti-) time-ordered products and (generalized) retarded products 
have in general no mutual coincidence regions in real energy-momentum space, which 
plagues the usual connection between the approach in real time and energy variables (RTF) 
(in terms of time-ordered products) and the use of the Green functions in complex energy 
variables. As a matter of fact, for the case of representations of field theory at finite 
temperature T = , the KMS-condition provides substitutes to the usual coincidence 
relations due to spectrum. These substitutes are relations which contain the Bose-Einstein 
factor (1 — e _/3w ) and thereby imply that coincidence relations can be only "essentially 
restored, up to exponential tails" in the limit of very high energies. 

Another consequence of the KMS-condition concerns the justification of the approach 
in purely imaginary time and energy variables. It is a standard practice of field-theorists to 
use the so-called "Euclidean n-point functions" , namely the Wightman functions at purely 
imaginary times (or Schwinger functions) and correspondingly the n-point Green functions 
at purely imaginary energies. In the vacuum representation of Quantum Field Theory, the 
justification of this practice is a by-product of the double analytic structure of the n-point 
functions which results from the interplay of spectral condition and locality. In particular, 
it is known that the Green functions at imaginary energies are obtained as the Fourier 
transforms of the corresponding Wightman functions at imaginary times. In the case of 
thermal representations of Quantum Field Theory, a similar double analytic structure of the 
n-point functions also results from the interplay of KMS-condition and locality. However, in 
a representation at temperature T = (3~ 1 , the n-point Wightman functions of the complex 
time-variables acquire periodicity conditions with period i(3 as a consequence of the KMS- 
condition. Correspondingly, the Green functions of the complex energy-variables (namely 
the Fourier-Laplace transforms of the generalized retarded functions) are expected to be 
completely determined from their values on a certain lattice of discrete purely imaginary 
energies of the form uj = ^p, with / integer, these values being the Fourier coefficients of 
the corresponding (z/3-periodic) Schwinger functions. 

These consequences of the KMS-condition have been understood to a large extent in 
the time-ordered path approach (see e.g. [4,7,9,23]). We shall exhibit them here in our 
general axiomatic setting and give complete results for the case of the two-point function, 
but only indicate how their generalization to the n-point functions could be worked out. 

The two-point function 

Substitutes to the coincidence relations in energy-momentum space: 

Let W(p) be the Fourier transform of the two-point function W(x). The KMS-condition 
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in the energy variable is expressed as a special case of formulas (5), (7) (see Sec.l iii')), 
namely: 

W(p) = e^W(-p) = Y^^J • (35) 

This way of writing Eq.(7) for the case of the two-point function (with the Bose-Einstein 
factor at the r.h.s. of this equation) is justified and unambiguous, as a relation between 
positive measures, in view of the additional postulate v) , namely, the time-clustering prop- 
erty for W. The latter implies that (after the addition of a suitable constant to the field) 
W(x) tends to zero when the time-variable xq tends to infinity (at fixed x) and corre- 
spondingly that the measure W(p) contains no term proportional to 5(u). Therefore, 
Eq.(35) determines uniquely the splitting of the spectral function C(p) = W(p) — VV(— p) 
and replaces the usual splitting ~W^ vac \p) = 9{uj)C^ vac \p) which results from the spectral 
condition in the vacuum representation. Eq(35) implies that C{p) and W(p) "coincide up 
to an exponential tail" at very high energies, while W(p) "vanishes up to an exponential 
tail" at negative energies of very high absolute value. 

The fact that C(p) does not vanish in general on any open subset of energy- momentum 
space (the examples given in Sec 2 are exceptions which will be commented in Sec. 6c 
below) implies that r(p) and a{p) do not coincide on any open set and therefore that 
the corresponding holomorphic functions in the tubes T + and T~ are not the analytic 
continuation of each other. 

Let f(p) be the Fourier transform of the time-ordered function t(x). In view of Eqs. 
(12), (13) and (35), the following relations hold: 

rip) = ~ia(p) + W(p) = KP [ _ e _ l /J , (36) 

which show that f and —if only "coincide up to an exponential tail" at very high energies. 
The double analytic structure: 

The KMS-condition (see Sec. 2 iii)) implies that there exists a holomorphic function 
W(zq,x) of zq in the strip — (3 < Ittlzq < whose boundary values on the edges of this 
strip are respectively: W{xq, x) = W(x) and W{xq — i/3, x) = W(x). Locality then implies 
that the holomorphic function W can be analytically continued as an z/3-periodic function 
of zq in the whole complex plane minus the cuts {zo = xq + iyo; \xq\ > \x\, yo = 1(3, 1 G Z}. 
Moreover, the jump AW of W across the right-hand cut xq > \x\ is (as in the case of the 
vacuum representation) : 

AW(x) = 9{xq) [W(x) - W(x)} = -ir(x). (37) 

If one now considers conjointly 
i) the Fourier-Laplace transform 

^o,p) = tAo/ dx f dx e^ koxo -^ r(x ,x) (38) 

y^r Jr s Jo 
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of r in its holomorphy domain Imko > 0, and 

ii) the Fourier coefficients of the (periodic) Schwinger function W(iy 0: x) : namely 



MP) = 77A2 I' d yo [ dxe-^Wiiyo^), (39) 

for I non-negative integer, one can show the following basic relations: 

Mp)=r(-^,P), 1 = 0, 1, 2,... (40) 

by means of a simple contour distortion argument. In fact, one checks that the integral of 
the product —(27v)~ 2 ie^ koZ °~P x ^ x W(z , x) over the complex cycle {z G [00 — z'f , —if] U 

[—if, if] U [if, 00 + if ]} x {x G 1R 3 } reduces to Eq (38) by shrinking the integration cycle 
onto the real set M. + x IR 3 , while it reduces to Eq (39) for the special values ko = ^p, / 
integer, / > due to the i/3-periodicity of W(z , x). 

Since f(ko,p) is of moderate growth in the complex half-plane Imko > (for each p), 
it can be characterized as the (unique) Carlsonian interpolation (see e.g. [34] p. 153) of 
the sequence of its values given by Eq.(40) at the set of discrete imaginary energies ^p, I 
integer, / > 0; a corresponding result holds for the advanced propagator a which satisfies 
in the lower half-plane equations similar to Eq.(40) with / < 0. This result together with 
formula (36) make completely clear the connection between the imaginary and real-time 
formalisms for the two-point function and their equivalence in the axiomatic framework. 

The n-point functions 

The derivation of various formulas which are substitutes to coincidence relations in energy- 
momentum space relies on the following fact. The relations (7) implied by the KMS- 
condition can be extended to a more general set of similar formulas, namely 

(1 - e-"*") (fy, ri(/)ri( = (fy, [ri(/), n(j)]fy), (4i) 

where 11(7) denotes the Fourier transform of any (product of) generalized retarded or (anti) 
time-ordered products of field operators depending on the energy-momentum variables 
Pi = (u>i,pi), % e /. Since all the distributions involved are tempered, it follows from 
Eq.(41) that all the expressions of the form (f2 ( g,n(J)n(J)f2 / g) decrease exponentially as 
e -/3|wj| £ or UJl tending to —00. 

By taking the latter property into account for all the terms of the sum at the r.h.s. of 
Eq.(20) (after applying a Fourier transformation to both sides of this equality), one easily 
checks the following statement: 

Consider for each cell a the region %, of (real) energy-momentum space which is defined 
by the condition u> = (cui, . . . ,u> n ) G 7 a ; then, in this region the corresponding distri- 
bution f a (p) "essentially coincides" with the Fourier transform f{p) of the time-ordered 
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expectation value {Cip, T({1, 2, . . . , n})^) up to terms which are exponentially decreasing 
with respect to the various energy variables inside 7 a . 

Formula (41) also applies to the r.h.s. of the following relations which result from 
Eq.(19) by a Fourier transformation: 



However, in contrast with what happens for the vacuum representation (in view of the 
spectral condition), formula (41) does not imply any support property for the commutator 
function at the r.h.s. of Eq.(42). Therefore, the holomorphic functions f ai (k) and r a2 (k) 
whose boundary values on the reals are linked by Eq.(42) do not admit mutual analytic 
continuation via the edge-of-the-wedge property as it is the case for the Green functions 
of the vacuum representation. 

The analytic structure of the n-point functions in the complex time-variables results 
from the KMS-conditions (see Eqs. (3), (4)) which one writes for all the permuted Wight- 
man distributions Wn(I, J), (I, J) denoting any pair of ordered sets forming an ordered 
partition of {1,2,..., n}. As in the case of the vacuum representation, there exists (for 

each n) a distribution Wn((zifi,xi),...,(z nj o,x n )) which is holomorphic in a union of 
permuted tubes in the space of complex time- variables z^o, 1 < i < n, together with 
regions of mutual crossover provided by locality. 

However, what is different and characteristic of the representations at temperature 
is the fact that the tubes have bounded polyhedral bases which form a multiperi- 
odic paving of the space of purely imaginary time-variables yo = (yi,o, ■ ■ • > Un,o)- In 
particular, the holomorphy domain of the contains all the Euclidean configura- 

tions Zj = (iyj^o,Xj), 1 < j < n, except those which satisfy conditions of the form 
Zj = Zk + ilfieo, I G Z, for all possible pairs (j,k). One therefore also expects a gen- 
eralization to take place concerning the double analytic structure in the time and energy 
variables which we have exhibited above for the two-point function. More precisely, one 
expects the derivation of formulas of the form (40) to be feasible in the general case; such 
formulas would relate the values of the generalized retarded n-point functions r a on ap- 
propriate lattices of purely imaginary energies to the Fourier coefficients of the Schwinger 
n-point functions which should be computed (for each r a ) on a corresponding periodicity 
pattern of yo-space. This is an algebraic problem which can be solved rather easily for 
n=3 and is expected to be solvable for general n with some amount of work [35]. 



a) Integral representation of the two-point correlation function 

It follows from Eqs. (14) and (35) that the Fourier transform of the two-point correlation 
function of the free field at temperature is 



f ai -r Q2 = (tip, [R a w,R ai j)]Ctp) . 



(42) 



6. COMPLEMENTS AND DISCUSSION 



1 e(u)5(u 2 -p 2 
2tt 1 - e~P w 




(43) 
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By multiplying both sides of Eq.(23) by (1 — e _/3w ) _1 and taking Eq.(35) into account, 
one then obtains the following general representation for the Fourier transform of the two- 
point correlation function of any local scalar field at temperature [3~ 1 : 

-I f POO 

W( u ,p) = -— T du dsp{u,s)W { /~ s \ujJ-u). (44) 

(27Tj2 JR3 Jo 



By coming back to the space-time variables, we obtain correspondingly the following 
integral representation for the correlation function itself, which is comparable to formula 
(21): 

/■OO 

W(x) = dm D(x;m) W { ™\x). (45) 
Jo 

Remark: 

The consideration of "generalized free thermal fields" completely defined by specifying their 
two-point function (and imposing the prescription that all the truncated n-point functions, 
with n > 2, should vanish) is now possible. We can say that all the information which is 
necessary to construct such a field satisfying all the general principles of Thermal Quantum 
Field Theory defined above (see Sec. 2) is encoded in the "weight-function" D(x;m) (or 
p(u,s)) of the integral representation (45) (or (44)) of its two-point function; the only 
restriction which the tempered distribution p(u, s) has to comply with is the following one: 
the convolution expressed (formally) by the double integral at the r.h.s. of Eq. (23) must 
give a positive measure. This is clearly the case if p(u, s) happens to be positive. 

b) Consequences of the relativistic KMS-condition 

The implications of the relativistic KMS-condition (presented above in Sec. 2, conditions 
iv), iv')) concern the existence of an analytic structure of the n-point functions in the 
complex space-time vector variables and the exploitation of corresponding exponential decay 
properties in the space of (real) energy-momentum variables. The latter are produced not 
only at very large negative energies but also at very large momenta; they in fact apply to 
energy-momentum regions and to distributions for which the relativistic spectral condition 
would imply vanishing properties (in the vacuum representation of field theory). 

For the case of the two-point function, a detailed study will be given in [32] with the 
following type of results. Even under the weakest form of relativistic KMS-condition (see 
[17]), the function W(zo,x) introduced in Sec. 5 admits an analytic continuation W(z) as 
a holomorphic function of the complex four- vector z = (zq, z) in a domain which contains 
the i/9-periodic "flat cut-domain" of W in (zq, x)-space. This shows (as a by-product) 
the regular (i.e. C°°) character of the distribution W(xq,x) with respect to the spatial 
coordinates x. Moreover, the distribution D(x,m) appearing as a "weight-function" in 
the representations (21) and (45) is also shown to have a C°°-dependence in x. Under the 
strongest form of relativistic KMS-condition (namely, the one presented as condition iv) in 
Sec. 2), W(z) is holomorphic in the union of the tube Tp together with all those which are 
obtained from the latter by the translations il(3eo, I G Z. All these tubes are connected 
together by complex neighbourhoods of the regions {z; yo = 1/3, x\ — x 1 < 0, / G Z} 
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given by locality. The weight-function D is also proved to admit an analytic continuation 
D(z, m) in a corresponding tube with imaginary basis (containing the origin) in the space 
of complex coordinates z. In the energy-momentum variables (oj,p), the product of the 

spectral function C(p) by the function e^ 1+ ^ )^~( 1 + aj ) 3 ] must be a tempered distribution 
and the "exponential tail" of f + if extends to the whole complement of the forward 
cone in IR 4 ; correspondingly, the "weight-function" p(u, s) of the representation (22) is 
exponentially decreasing with respect to \u\. Similar properties could be derived for the 
n-point functions. 

c) Feynman-type operations on two-point functions and discussion of examples 

The double analytic structure of thermal two-point functions which we have described (see 
Sec. 5) can be shown to be preserved under the two basic Feynman-type operations, namely 

i) "N-line Wick-contraction" (represented by the diagram with two vertices z\ and z 2 con- 
nected by N lines): 

With such a diagram V is associated the product W^[r](^i — z 2 ) of N two-point functions 
taken at the same complex point z\ — z 2 (each factor W^(z\ — z 2 ) being associated 
with a line A). It is clear that such a product of holomorphic functions still satisfies 
the same properties (z/3-periodicity, analyticity domain expressing the relativistic KMS- 
condition, boundary values on the reals in the sense of distributions, locality) as each 
individual factor. These properties characterize Wp] as being a thermal two-point function 
whose associated propagator f^(k) or a[p](/c) can be computed on the discrete sequence 
of energies ui = , / G Z, (see Eq.(40)) as a (discrete) convolution of the N propagators 
f(>) or d(>) taken on the same sequence of imaginary energies. Moreover, the Fourier 
transform of the correlation function VV[r] (p) is also the convolution product on the (real) 
energy-momentum space of the Fourier transforms of the N correlation functions VV(a) • 

ii) "Vertex convolution": 

The vertex convolution of the thermal two-point functions and W( 2 ) is defined as the 
following convolution product on Euclidean space-time: W(z\ — z 2 ) = (W(i) * W(2)) (z± — 

z 2 ) = J 2 (3 dzo J R3 dz W(i)(zi — z) W(2){z — z 2 ). The resulting holomorphic function W 

still satisfies all the structural properties of thermal two-point functions and the associ- 
ated propagator f(k) in complex energy-momentum space is the ordinary product of the 
propagators f\{k) and f 2 (k) associated with W\ and W 2 . 

In particular, if one starts from the free thermal two-point functions (see Sec.3, Ex- 
ample 1, and Sec. 6a)), one can construct as an application of the operation i) appropriate 
perturbative two-point functions associated with "self-energy bubble diagrams" of the 
thermal interacting field. Then, the operation ii) defines in a rigorous way the two-point 
function obtained by "resummation of the corresponding iterated self-energy diagrams". 
If the self-energy contribution is denoted (in complex energy-momentum space) by S(/c), 
the associated "complete" propagator is given, as usual, in the complex domains (T+ and 
T") by the formula A(k) = ^iZ^Z^y 
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All these structural properties could of course be derived similarly for non-scalar fields 
such as those of QCD. One can thus understand in this general framework such computa- 
tions of gluon propagators as those given by Weldon and Pisarski [5,6] (see Sec. 3, Example 
2). In the latter, however, a high-temperature approximation has been taken which makes 
the expressions (16) and (17) of the propagators A; and A t somewhat peculiar, as far as 
the support and decrease properties of the associated spectral functions are concerned. 

We first notice that, apart from those of the free field and of the two-lined bubble 
diagram, the spectral functions of all the perturbative two-point functions generated by 
the previous operations i) and ii) have no support restrictions; however, they all enjoy 
a property of exponential decrease at large momenta which expresses the fact that the 
relativistic KMS-condition is satisfied. These two properties are violated by the spectral 
functions of the WP-propagators, since: 

a) the latter only satisfy a condition of power decrease at large momenta: but this "hard 
thermal loop approximation" is generally used only in the low momentum region; 

b) in the usual relativistic spectral region \u>\ > \p\, their support is restricted to a pair of 
hypersurfaces of the form u = ±ut,i(p) which correspond to sharp terms containing the 
factors S(u> =F ^t,i(p))', the existence of such real dispersion laws, interpreted as "plasmon 
modes" , calls for further comments which open our last topic. 

We conclude these remarks by noting that a general perturbative approach to the con- 
struction of thermal correlation functions, based on the systematic exploitation of locality 
and KMS-condition, has been proposed by Steinmann [36]. 

d) Representations of particles 

The concept of "dispersion law" u> = uj part {p) for representing a "particle" (or "mode") 
in Thermal Quantum Field Theory may seem the most natural one to be inherited from 
the familiar formalism of quantum field theory in the vacuum representation, having in 
mind that the hypersurface defined by this law in the space of the energy-momentum 
variables might not be (in general) a relativistic hyperboloid shell. Such a hypersurface 
should appear as the singular set of a real pole in the thermal propagators of the theory, 
or equivalently as the support of a 5 -term in the corresponding spectral functions. How- 
ever, this attractive picture turns out to be false, except precisely in the trivial case of 
a relativistic free particle moving across the thermal bath without any interaction: so is 
the content of the Narnhofer-Requardt-Thirring Theorem [37] which has been proved by 
these authors in the general framework of "Local Quantum Physics" [22] . This frustrating 
result holds true, even if one allows the support of the 5-term to be imbedded in a region 
where the spectral function has a non-zero continuous background. As a matter of fact, 
this is not so surprising since a 5-term in C(p) would imply that the correlation function 
W(t, x) has the "normal" slow decay property (as |t|~2) along any world-line x = vt, and 
therefore that the particle is not submitted to dissipative effects due to the interactions 
with the thermal bath. For example, the WP-gluon-propagators considered above repre- 
sent an approximation which does not provide a realistic description in terms of particles 
in a thermal equilibrium state. 
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A reasonable way out of this deadlock, which has been proposed in particular by 
Landsman [38] , consists in adopting the same viewpoint as for the representation of unsta- 
ble particles in the vacuum state of field theory. This amounts to assume that the retarded 
propagator (or a distinguished part of it) admits an analytic continuation from the upper 
/co-plane (or the tube T + ) into the lower /c -plane, across (part of) the reals, and that a 
complex pole ko = uJ pa rt(p) ~ ilpart{p) is present in this second-sheet domain. Of course, 
this complex dispersion law should present some characteristics which would distinguish 
the type of decay (due to statistical dissipative effects) of a particle in a thermal bath from 
the one (due to intrinsic unstability) of a resonance. In particular, one would expect the 
law to be such that the dissipative effects are barely felt by the particle at rest, namely 
that W(t, vt) keeps the \t\~^ behaviour for v = 0, while being exponentially decreasing in 
t for v ^ 0. This is certainly not the case for the simplest ways of choosing the "width" 
Ipartip) (i-G- constant or proportional to ujpart 

At this point, we would like to advocate (as in [18,19]) that our general integral 
representations (23) and (45) of the thermal spectral functions and correlation functions 
provide a somewhat more natural prescription for representing particles: let us assume 
that, as in the Kallen-Lehmann formula for the spectral function in the vacuum state, 
a particle is associated with a 6 -term in the weight-function p(u, s) of the representation 
(21), namely, with a term of the form p par t{u) S(s — Too)- The representation (45) then 

contains the distinguished term D part (x) Wp°\x), where D part (x) is (up to a factor) the 
Fourier transform of p par t(u). It is now clear that, if we interpret D part (x) as a "dissipative 
(or damping) factor" which should decrease, for instance, exponentially at large x, this 
distinguished term is a good candidate for representing a particle behaviour. (Note that 
it would satisfy the exponential decrease property along all world-lines x = vt, but would 

3 ( \ 

behave as \t\~z at rest, like the free correlation function W^ .) Since in this case, the 
function p par t would be holomorphic in a tube domain of the complexified variables u, one 
can prove that the corresponding propagator represented by formula (34) would then have 
analytic continuation properties in the lower half-plane of the variable ko- This indicates 
(and a further general study confirms it as well as the special example presented in [18,19]) 
that this 5 -representation of particles, far from being contradictory with the previous one 
based on "second sheet complex poles" of the propagators, selects inside the latter class 
candidates which might be the most appropriate representatives of the notion of particle 
in Relativistic Thermal Quantum Field Theory. 

Finally, one can also show that the integral representations (21) and (23) are a useful 
tool for investigating the manifestations of "thermal Goldstone particles" in the case of 
"spontaneous symmetry breaking" ; this will be the object of a further work [39] . 

In conclusion, we have presented a general framework "a la Wightman" for the con- 
cepts of Thermal Quantum Field Theory and a set of (preliminary) results, mainly ex- 
pressed in terms of analytic structural properties of the n-point functions of the fields. 
These results are already rather complete for the case n = 2. They are useful for clarifying 
some structural aspects of perturbative results, obtained in the previous years by various 
authors. Moreover, the general integral representation of thermal two-point functions that 
we have obtained is interesting under several respects: providing a complete knowledge 
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of "generalized free thermal fields", it also suggests a rather promising approach to the 
concept of particle, integrating the previously known ideas and consistent with the general 
principles of relativistic thermal field theory. Many open problems remain, as far as the 
properties and the use of n-point functions are concerned (for n > 2); among the important 
ones, let us mention in particular the problem of the possible definition of several-particle 
states as a manifestation of the field interactions in a spirit comparable to the collision 
theory in the vacuum state. 
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